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Abstract 

We consider equations of the form 

(d^u) (t, x) - Au(t, x) = f(t, x), t>0,x£ R n , 

where is a distributed order derivative, that is 

1 

o 

B*") is the Caputo-Dzhrbashyan fractional derivative of order a, fi is a positive weight 
function. 

The above equation is used in physical literature for modeling diffusion with a loga- 
rithmic growth of the mean square displacement. In this work we develop a mathematical 
theory of such equations, study the derivatives and integrals of distributed order. 



Key words: distributed order derivative; distributed order integral; ultraslow diffusion; 
fundamental solution of the Cauchy problem. 
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1 INTRODUCTION 



Fractional diffusion equations with the Caputo-Dzhrbashyan fractional time derivatives 

(d! q) w) (t,x) -Bu(t,x) = f(t,x), i>0,i6i n , (1.1) 

where < a < 1, B is an elliptic differential operator in the spatial variables, are widely used 
in physics to model anomalous diffusion in fractal media. 

Physically, the most important characteristic of diffusion is the mean square displacement 



(Ax) 2 = J \x-Z\ 2 Z(t,x-t)d{i 

of a diffusive particle, where Z is a fundamental solution of the Cauchy problem for the diffusion 
equation. In normal diffusion (described by the heat equation or more general parabolic equa- 
tions) the mean square displacement of a diffusive particle behaves like const ■ t for t — > oo. A 
typical behavior for anomalous diffusion on some amorphous semiconductors, strongly porous 
materials etc is const • t a , and this was the reason to invoke the equation (1.1), usually with 
B = A, where this anomalous behavior is an easy mathematical fact. There are hundreds of 
physical papers involving equations (1.1); see the surveys [231 123] - The mathematical theory 
was initiated independently by Schneider and Wyss [37] and the author [THEE]; f° r more recent 
developments see [91 [101 [15] and references therein. 

A number of recent publications by physicists (see 0, HI EJ [25l EH] and references there) is de- 
voted to the case where the mean square displacement has a logarithmic growth. This ultraslow 
diffusion (also called "a strong anomaly" ) is encountered in polymer physics (a polyampholyte 
hooked around an obstacle), as well as in models of a particle's motion in a quenched random 
force field, iterated map models etc. In order to describe ultraslow diffusion, it is proposed to 
use evolution equations 

(©^w) (t,x) -Bu{t,x) = f{t,x), t>0,xeR n , (1.2) 

where is the distributed order derivative (introduced by Caputo [2J) of the form 

l 

(B^V)(0 = y"(© (a V)WM«)da, (1.3) 
o 

fi is a positive weight function. 

The above physical papers contain some model calculations for such evolution equations. 
There are only two mathematical papers on this subject. Meerschaert and Schefffer [22] devel- 
oped a stochastic model based on random walks with a random waiting time between jumps. 
Scaling limits of these random walks are subordinated random processes whose density func- 
tions solve the ultraslow diffusion equation. The solutions in [22] are understood as solutions 
of "algebraic" equations obtained if the Laplace transform in t and the Fourier transform in x 
are applied. 



3 



Umarov and Gorenflo [39] applied to equations (1.2) Dubinskij's theory [8] of analytic 
pseudo-differential operators. This leads to solvability results for (1.2) in the spaces of an- 
alytic functions and dual spaces of analytic functionals. Such a theory is very different from 
the theory of parabolic equations; results obtained this way "do not feel" the difference between 
the equations (1.2) with B = A and B = —A. 

The aim of this paper is to develop a theory of the model equation (1.2) with B = A 
comparable with the classical theory of the Cauchy problem for the heat equation. In particular, 
we construct and investigate in detail a fundamental solution of the Cauchy problem for the 
homogeneous equation (/ = 0) and the corresponding kernel appearing in the volume potential 
solving the inhomogeneous equation, prove their positivity and subordination properties. This 
leads to a rigorous understanding of a solution of the Cauchy problem - it is important to know, 
in which sense a solution satisfies the equation. In its turn, this requires a deeper understanding 
of the distributed order derivative (1.3), the availability of its various forms resembling the 
classical fractional calculus [36]. We also introduce and study a kind of a distributed order 
fractional integral corresponding to the derivative (1.3). A Marchaud-type representation of 
the distributed order derivative (based on a recent result by Samko and Cardoso |35J) is the 
main tool for obtaining, in the spirit of [191 E] , uniqueness theorems for the Cauchy problem 
for general equations (1.2) in the class of bounded functions and, for n — 1 and B = d 2 /dx 2 , 
in the class of functions of sub-exponential growth. 

Comparing with the theory of fractional diffusion equation (1.1) we see that the distributed 
order equations (under reasonable assumptions regarding /i) constitute the limiting case equa- 
tions, as a — > 0. That is readily observed from estimates of fundamental solutions having, 
as |sc| — ► oo, the estimate exp ^— a|x|2^j, a > 0, for the fractional diffusion equations, and 

exp(— a\x\) in the case of ultraslow diffusion. 

In fact, we begin with the "ordinary" equation O^w = Am, A G R. If A < 0, already this 
equation demonstrates a logarithmic decay of solution at infinity; see Theorem 2.3 below. 

In general, the theory presented here is an interesting example of subtle analysis (with 
kernels from L\ belonging to no L p , p > 1, etc) appearing in problems of a direct physical 
significance. 

2 Distributed Order Derivative 

2.1. Definitions. Recall that the regularized fractional derivative of a function (p G C[0, T] 
(also called the Caputo or Caputo-Dzhrbashyan derivative) of an order a G (0, 1) is defined as 




< t < T, 



(2.1) 



o 



if the derivative in (2.1) exists. If <p is absolutely continuous on [0,T], then 



t 




(2.2) 



o 



(see [9]). 
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Let //(£), < t < 1, be a continuous non-negative function, different from zero on a set of 
a positive measure. If a function if is absolutely continuous on [0, T], then by (1.3) and (2.2) 



t 



(© (/x V) (t) = I k(t- t)<p'(t) dr (2.3) 

where 



o 



l 

k(s) = I— -fj,(a)da, s > 0. (2.4) 



ri-Q 







It is obvious that k is a positive decreasing function. 
Note that for an absolutely continuous ip, 



t 

d 



J k(t- r)(p(r) dr= j k{s)^{t - s) ds = J k(s)<p'(t - s) ds + k(t)<p(0), 



so that 



(D ( "V) (*) = j I k{t - t)<p(t) dr - k(t)(p(0). (2.5) 
o 

The right-hand side of (2.5) makes sense for a continuous function (p, for which the derivative 
d 1 

— j kit — t)(p(t) dr exists. Below we use (2.5) as a general definition of the distributed order 
derivative D" 1 ^. 

The necessity to use the regularized fractional derivatives, not the Riemann-Liouville ones 
(defined as in (2.1), but without subtracting t~ a <p(0)), in the relaxation and diffusion problems, 
is caused by the fact that a solution of an equation with a Riemann-Liouville derivative typically 
has a singularity at the origin t = (see, for example, [H]), so that the initial state of a 
system to be described by the equation is not defined and requires a regularization. However, 
mathematically such problems are legitimate. A distributed order derivative with a constant 
weight, based on the Riemann-Liouville fractional derivative, was introduced by Nakhushev 
(see (26J). The diffusion equation with such a time derivative and a single spatial variable 
was investigated by Pskhu [5U]. For other definitions of variable order and distributed order 
derivatives see also [T7J [21] and references therein. In this paper we use only the derivatives 
(2.1) and (2.5). 

2.2. Asymptotic properties. Since the kernel (2.4) is among the main objects of the 
distributed order calculus, it is important to investigate its properties. 

Proposition 2.1. If (i e C 3 [0, 1], ^ 0, then 

k(s) ~ s" 1 (logs)- 2 /i(l), s^O, (2.6) 
k'(s) ~ -s~ 2 (logs)~ 2 /i(l), s^O. (2.7) 
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Proof. Denote r = — logs (— > oo), ip(a) = - . Then 

T(l — a) 

1 

k(s) = J tP(a)e ra da. 
o 

Integrating twice by parts, we get 

i 

jfe(s) = r - 2 / ^"(a)e ra rfa - r"V(0) - r~ 2 [^'(l)e r - V'(0)] . 



We have 

,., , «'(a)r(l — of) + //(a)r'(l — a) 

^ (a) = pa^f ' 

so that ^'(l) = — /i(l), and another integration by parts yields the relation 

k(s) = //(l)r~V + 0(r~V), r -> oo, 

which implies (2.6). The proof of (2.7) is similar. ■ 

It follows from (2.6) that k G Li(0, T); however k ^ Lp for any /3 > 1. Note also that one 
cannot integrate by parts in (2.3) because, by (2.7), k! has a non-integrable singularity. 
Throughout this paper we use the Laplace transform 



fC{p) = J k(s)e- ps ds, Rep > 0. 
o 

Using (2.4) and the relation 

oo 

s e p ds = 



pl-a 



(see 2.3.3.1 in [28]), we find that 

i 

/C(p) = J p a -V(a)rfa. (2.8) 
o 

It will often be useful to write (2.8) as 

i 

K{p) = [ e^ 1)logp fx(a)da. (2.9) 



Taking the principal value of the logarithm we extend JC(p) to an analytic function on the whole 
complex plane cut along the half- axis R_ = {Imp = 0, Rep < 0}. 
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Proposition 2.2. (i) Let fi G C 2 [0, 1] . J/peC\R_, \p\ -> oo, t/ien 



^) = ^ + 0((logH)- 2 ). (2.10) 

More precisely, if /i G C 3 [0, 1], t/ien 

/C(P) = - + O ((log b|)- 3 ) . (2.10') 

iogp (logpj z 

(ii) Lei // G C[0, 1], /x(0) ^ 0. Ifp G C \ R_ ; p -> 0, then 

K(p) ~ p- 1 flogi) /i(0). (2.11) 

(iii) Let // G C[0, 1], /i(ct) ~ aa A ; a > 0, A > 0. Ifp eC\t_,p^0, then 

}C(p)^aT(l + X)p- 1 (\og-) . (2.11') 



Proof, (i) Integrating by parts, as in the proof of Proposition 2.1, we find that 



i 



y e a >(a) da = + O (r' 2 e r ) , r -> oo, 



o 



which implies (2.10). The relation (2.10') is proved similarly. 

(ii), (iii). The relations (2.11) and (2.11') follow from the complex version of Watson's 
lemma ([27], Chapter 4). ■ 



In some cases it is convenient to use a rough estimate 

|/CG9)|<C1jWlog^) , \p\< Po , 



IpI 

valid for any /i G C[0, 1]. This estimate follows from general results about the behavior of the 
Laplace transform near the origin (see Chapter II, §1 of [7j). 

2.3. "Ordinary" equations. Let us consider the simplest equation with a distributed 
order derivative, that is 

(D M «J (t) = \u x {t), t > 0, (2.12) 

where A G R, and it is assumed that a solution satisfies the initial condition u(0) = 1. A 
solution of (2.12) should be seen as an analog of the exponential function t i— > e xt of the classical 
analysis and the function t i— > E a (Xt a ), where L7 Q is the Mittag-Leffler function, appearing for 
the equation with the regularized fractional derivative of order a G (0, 1) (see [9]). The equation 
(2.12) with A < is discussed in [16J as the one describing distributed order relaxation. The 
uniqueness of a solution will follow from the uniqueness theorem for the equation (1.2); see 
Theorem 6.1. Of course the method of proof of the latter theorem can be used to prove 
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separately the uniqueness for a much simpler equation (2.12). Below we assume that /i G 
C 2 [0, 1], /i(l) ^ 0, A ^ 0; evidently, u (t) = 1. 

Applying formally the Laplace transform to the equation (2.12) and taking into account the 
initial condition u(0) = 1, for the transformed solution u\{p) we get 

The function (2.13) is analytic on the half-plane Rep > 7, if 7 > is large enough. We have 
u\{p) ~ p~\ p = cr + ir, a, t G R, |r| ^00. Therefore [7] tZX is indeed the Laplace transform 
of some function U\(t), and for almost all values of t, 

7+ioo 

UA (t) = _ / \> dp. (2.14) 

AV ; di27Ti jf p p/C(p) - A v ; 

7— ioo 

Let I < a> < 1. We will often use the contour in C consisting of the arc 

T^u; = {p G C : |p| = 7, 1 argp| < ^vr}, 

and two rays 

r 7,^ = {P e C : I argp| = ^vr, |p| > 7}, 

= {P e C : I argp| = -uir, |p| > 7}, 

The contour S 1}U1 is oriented in the direction of growth of argp. 
By Jordan's lemma, 

d 1 /" e pt Up) 



dtl-ni J p p/C(p) — A 
In contrast to (2.14), here we may differentiate under the integral, so that 

1 / ,„< K M 



-*w=2sy c 'sfa*' (2 - i5) 

Note that 7 is chosen in such a way that p/C(p) 7^ A for all p G (this is possible, since 
p/C(p) — >■ 00, as |p| — > 00). 

The next result establishes qualitative properties of the function u\ resembling those of the 
exponential function and the Mittag-LefHer function. Recall that a function u G C°°(0,oo) is 
called completely monotone [13], if (— l) n u^ n \t) > for alH > 0, n = 0, 1,2,.. .. 

Theorem 2.3. (i) The function u\(t) is continuous at the origin t and belongs to C°°(0, 00). 
(ii) If A > 0, then U\(t) is non- decreasing, and u\(t) > 1 for all t G (0, 00). 
(hi) If X < 0, then u\(t) is completely monotone. 
(iv) Let A < 0. Iffi(0) + 0, then 

u x (t) ~ C(logt)~\ t^oo. (2.16) 

If yu(a) ~ aa u , a — > (a > 0, v > 0), then 

u\(t) ~ C(logt)- 1 -^, t^oo. (2.17) 

Here and below C denotes various positive constants. 
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Proof. The smoothness of u\ for t > is evident from (2.15). The integral in (2.15) is the 
sum of the integral over T ltLU (obviously continuous at t = 0) and the functions 



A w 2m J pK{p) — A 



We have 



«a (*) + «a (*) = - Im S / e' re . „ v , . ; r dr 

7T I J re^ 7r /C(re^ 7r ) - A 



— Im / r _1 e* rel dr H — Im / — - — : — — — : — r — dr. 



7r J 7r J r (re iuJn JC(re iuJ7T ) - A) 

7 7 

The second summand is obviously continuous at t = 0. The first summand equals 

oo oo 

- / r~ 1 e trcosujn sm(trsmL}iT)dr = - / s _1 e~ s sin(-stanu;7r) ds 

7 — 7tCOSOJ7T 

(recall that cos urn < 0), and this expression is continuous at t — 0. 

Let A > 0. The function p i— > is the Laplace transform of the function £ i— > e A *. 

p — A 

Therefore [13] it is completely monotone. The Laplace transform of the function u' x (t) equals 

—/ v / n x pHp) , A 



pfC(p) — A pfc(p) — A 

(strictly speaking, we have to use this formula, together with the asymptotics of JC{p), to prove 
the existence of the Laplace transform of u' x ). 
On the other hand, the function 

i 

pK{jp) = / p a fj,(a) da 



is positive, while its derivative is completely monotone. By the Criterion 2 of the complete 
monotonicity (see Chapter XIII of [13]), the Laplace transform of the function u' x {t) is com- 
pletely monotone. It follows from the Bernstein theorem about completely monotone functions 
and the uniqueness property for the Laplace transform that u' x (t) > for all t > 0, whence u x 
is non-decreasing and u x {t) > 1. 

Let A < 0. Up to now, 7 was chosen so big that p)C(p) 7^ A for all p G S 1;UJ . In fact, 

1 

lmplC(p) = / \p\ a fj,(a) sin(a argp) da, 



so that Imp/C(p) = only for axgp = 0. Meanwhile, if argp = and A < 0, then pK,(p) — A = 
i 

/ \p\ a fi(a) da — A > 0. Therefore, the above integral representation of u\ holds for any 7 > 0. 


Since \JC(p)\ < C|p| _1 (log A) for small \p\, we find that 



K{p) 



pK{p) - A 



< ( '\l>\ (log |^| 



whence 



pK{p) -\ dp 



Cei 1 
< .-0, 



logi 



as 7 — > 0. 

Considering other summands in the integral representation of i^, we see that 



— 1 tre %wl ' 

- I r V re 



dr — 

7T 



s 1 e s sin(stanu;7r) ds 



-— s e s sin(s tanc<J7r) ds, 

7T 



7 —7* COS 

as 7 — > 0. Next, we have to consider the expression 



00 

— / Im 

7T 



Re 



1 



re iuJ7T IC(re iuJ7T ) - A 

00 

+ ~J Re 

7 



Im 



re ia;7r /C(re ia;7r ) - A 



dr^h + h. 



We have 



n tre 1 ' 



Im 



j,— lgtr cos am 



sin(tr sincj7r), 



and this expression has a finite limit, as r — > 0. Since also pJC(p) — > 0, as p — > 0, we see that 
we may pass to the limit in 7i, as 7 — > 0. 

In order to consider J 2 , we have to study the function 

$(r, a;) = Im ■ 



re^ 7r /C(re^ 7r ) - A' 



a(s— icon) 



We have 

1 



s = — logr — > 00, as r — > 0, so that 

1 

f e~ as sin(aw7r)/i(a:) da 

$(r, a;) = - 



/z(a) da, 



f e as cos(aa;7r)/i(a) da — A 



+ 



-1 2 • 



J e as sin(ac<J7r)/z(a) da 
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As s — > oo, the denominator tends to A 2 , while in the numerator 
1 1 



e as sm(auir)u(a) da < C / ae as < — — , 

(logr)^ 



C 



c 



o o 
This makes it possible to pass to the limit in I2, as 7 — > 0, so that 



u x {t) 



s- 1 e- s sm(staiiLj7r)ds + - I r ' l e trcos ^ 



A 



_| I r -l g ircosw7r 

7T 



sin(tr sin u;7r)\I/(r, u) dr 

cos(tr sin cJ7r)$(r, a;) dr (2.18) 



where 



f e as cos(aco>7r)/z(a) da — A 



J e~ as cos(au;7r) /x(a) da — A 



+ 



n 2 ' 



J e~ as sin(aa;7r)/i(a) da 



s = — logr. 

In (2.18), we may pass to the limit, as u — ► 1. It is easy to see that the first two terms in 
(2.18) tend to zero, so that 



u x (t) = ^ J r-V* r $(r, 1) dr, 



(2.19) 



J r a sm(air) n(a) da 



d>(r, I) 



J r a cos(a7i)fi(a) da — A 



+ 



J r a sin(a7r)/x(a) da 



2 • 



Since A < 0, it is seen from (2.19) that is the Laplace transform of a positive function. 
Therefore u\ is completely monotone. 

Let A < and /x(0) 7^ 0. As we have proved, u\ is monotone decreasing. It follows from 

C 

(2.13) and (2.11) that u\(p) ~ r, p — > +0. Applying the Karamata-Feller Tauberian 

plog- 

theorem (see Chapter XIII in [13J) we get (2.16). Similarly, if //(a) ~ aa u , a — > 0, we use the 
asymptotic relation (2.11'), and the same Tauberian theorem yields (2.17). ■ 

A non-rigorous "physicist-style" proof of the statement (iii) was given in [16J, where the 
asymptotics (2.16) was also found for the case /i(a) = 1. 



3 Distributed Order Integral 

3.1. Definition and properties. Suppose that W^'u = f, u(0) = 0. Applying formally the 

Laplace transform u >— > u, we find that u(p) = \ f(p)- The asymptotic properties of Kip) 

pJC(p) 

11 



show (see [7]) that the function p \— > is the Laplace transform of some function x(t), and 

plC(p) 

7+ioo 

d 1 f e pt 1 



* {t) = 7t^l / 7'i«W* 7> °- (31) 

7— ioo 

It is natural to define the distributed order integral 1^, as the convolution operator 

t 

{l M f)(t) = [ x(t-s)f(s)ds. 



Proposition 3.1. Suppose that fi e C 3 [0,1] ; 7^ ; and either fi(0) / or M a ) ~ aa v , 
a > 0, v > 0. Then: 

(i) x e C°°(0, 00) ; and x is completely monotone; 

(ii) /or small values of t, 

x(*)<Clogi (3.2) 

1x^)1 <cr x log i (3.3) 

Proof. As in Sect. 2, we deform the contour of integration in (3.1) and differentiate: 

"® = kj*m* (3 - 4) 

We will need information about the asymptotic behavior of . By (2.10'), 

/C(p) = ^-/^ + c(p), c(p) = O f — -L^) , |p|-oo. 
logp (logp) 2 \( lo g|Pl) / 

Then we can write 

1 logp M(l) -/i(l)c(p)(logp) 3 + [^(l)] 2 - M(l)c(p)(k>gp) 2 

/c(p) Mi) Ki)] 2 [Mi)] 2 [/i(i)iogp-M(i) + c(p)(io gP ) 2 ] ' 

whence 

1 logp ;/(!) , n / 1 \ 



/c(M Mi) [Mi)] 2 V(iogbl);' 

The integral in (3.4) consists of the integral over T lyW (a function from C°°[0,oo)) and 

<± 

7,i 



integrals over . Each of the latter ones is estimated, due to (3.5), by an expression 



C I er art r- x log r dr ~ Clog -, t -> 
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(a,C > 0; see the asymptotic formula (13.49) in [32]). This implies (3.2). The proof of (3.3) is 
similarly based on the same asymptotic relation from [32] . 

In order to prove that x is completely monotone, we proceed as in the proof of Theorem 
2.3, to transform (3.4) into a representation by a Laplace integral. First we pass to the limit, 
as 7 — > 0. This is possible because, by Proposition 2.2, either 



1 ~ ^(O)log-, p^O, 



if /i(0) ^ 0, or 



pK(p) C ( l0g p 



1+v 



, p^O, 



(3.6) 



(3.7) 



if h(ol) ~ aa u , a — > 0. Both the relations (3.6) and (3.7) are sufficient to prove that the integral 
over T 7jW tends to 0, as 7 — 
come to the representation 



over T 7jW tends to 0, as 7 — * 0, while the 7 — > limits of both the integrals over r±, exist. We 



1 



Im < e iW7r / e tr£l ' 



TV 



dr 



re iu)7T IC(re iu)7T ) 



(3.8) 



We find, introducing the parameter s = — logr — > 00, as r — > 0, that 



1 1 




e as (cos(acj7r) + z sin(ac<j7r)) / u(a) da. 



Taking into account the logarithmic behavior of the integrand of (3.8) near the origin, we may 
pass to the limit in (3.8), as uj — > 1, and we get that 



00 

Kit) = l -je 



■tr 



J r a sm(a7r)fi(a) da 




n 2 



J r a cos(aTr)fi(a) da 



+ 



J r a sin(a7r)/^(a) da 



2 



as desired. ■ 

Note that, by (3.2), x G 1^(0, 00). 

3.2. The Marchaud form of the distributed order derivative. If / G Li(0,T), 



1^/, then u = »c*f, 



(t) = * x * /)(*) = ^(1 * /) = ^ / f(r) dr = f(t) 



dt 



dt 
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almost everywhere. Thus IJWIW = / on Li(0, T). 

The identity (k * x)(t) = 1 (almost everywhere), which follows from the fact that the 
product of the Laplace transforms tC(p) and pf ^^ equals -, means that x is a Sonine kernel (see 
[35J). Since both the functions k and x are monotone decreasing (obviously, k is completely 
monotone), we are within the conditions [22], under which the operator D^', on functions 
u = I^f, f G L p (0, T), 1 < p < oo, can be represented in the form 

(B^u) (t) = k{t)u{t) + J k'{r)[u{t - t) - u{t)) dr, < t < T, (3.9) 

o 

where the representation (3.9) is understood as follows. Let 

(*e«) (t) - 



f k'(r) [u{t - r) - u{t)} dr, if t > s, 
0, if < t < e. 



Then 

lim || (D^u) (t) - k(t)u(t) - (* e u) W|| Lp(0iT) = 0. (3.10) 

The representation (3.9), similar to the Marchaud form of a fractional derivative [32], will 
be useful for our proofs of uniqueness theorems, because the integral operator in (3.9) has the 
form enabling the maximum principle approach. On the other hand, the precaution we made 
understanding (3.9) in terms of (3.10) cannot be easily avoided due to a strong singularity of 
k' in accordance with the asymptotics (2.7). 

4 The Equation of Ultraslow Diffusion 

4.1. A fundamental solution of the Cauchy problem. Let us consider the equation (1.2) 
with B = A, that is 

(t, x) = Au{t, x), xeR n ,t>0. (4.1) 

In this section we construct fundamental solution Z(t,x) of the Cauchy problem, a solution of 
(4.1) with Z(0, x) = 6(x), and obtain its estimates. 

Below we use the following normalization of the Fourier transform: 

u(£) = J e tx '^u(x) dx, 

i" 

so that 

For a radial function u(r), r = \x\, 

oo 

u(r) = 2W 2 (0 1_f f p^u{p)J^ 1 {rp)dp 1 (4.2) 



o 
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where J v is the Bessel function. 

Applying formally the Laplace transform in t and the Fourier transform in x, we find that 



By (4.2), 



oo 

n/2 



Z(p,x) = (2tt)"§|x| 1 -5/C(p) J - p7 - s -—^jn_ 1 (\ x \s)ds. (4.3) 







plC(p) + s s 



It is known ([29J, 2.12.4.28) that 



CO 



U 1 ;Ucy) dy = z v K u {cz), -I < v <%, (4.4) 







y 2 + z 2 2 
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where If^ is the McDonald function. If n < 4, then the above restriction upon v — 1 is 

satisfied, and (4.4) implies the representation 



Z(p,x) = (27r)-f |x| 1 -f/C(p)(p/C(p))2(f- 1 )K f _ 1 (|x| V ^(p)). (4.5) 
We have simpler formulas in the lowest dimensions - if n = 2, then 

Z(p,x) = ^-JC(p)K {\x\^/pX(p)y, (4.6) 

if n = 1, then 

Zfa X ) = lJ^=e~^V^) (4.7) 



because K-x/ 2 (z) = K 1 / 2 (z) = z (see pQ). 

The function K v decays exponentially at infinity: 




K v {z) ~ A /— e \ z^> oo, 



while K v {z) ~ as 2 — > (if f > 0), and K (z) ~ — log 2;. We see that the function on 

the right in (4.5) belongs to LiiW 1 ) in x for any n, not only for n < 4. Using the identity 



rJ v (hr)K v {cr) dr = b v c~ u {b 2 + c 2 )" 1 







([29], 2.16.21.1) we check that the inverse Fourier transform of the right-hand side of (4.5) 
coincides with K.(p) (pK.(p) + |£| 2 ) . Therefore the formula (4.5) is valid for any n. 

Let us consider estimates of the function Z and its derivatives. Qualitatively, the behavior 
of Z is similar to that of the fundamental solution of the Cauchy problem for the fractional 
diffusion equation (1.1) (see [19], [TO] E] ) . In addition to the singularity at t — 0, Z(t,x) has, if 
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n > 1, a singularity at x = (a logarithmic singularity, if n — 2, and a power one, if n > 3). 
As usual Z(t,x) — > as £ — > 0. This means that the singularity at t — becomes "visible" 
near the origin in x. In fact, we obtain separate estimates for a small \x\, showing the character 
of singularities in t and x, and for a large In addition, subsequent applications of the 
fundamental solutions require estimates of O^Z, applicable simultaneously for all i^O, and 
uniform in t. Of course, estimates for O^Z at the origin and infinity can be obtained from the 
relation D^Z = AZ. 

All the above estimates deal with a finite time interval, t G (0, T], and it is this kind of 
estimates, that is needed to study the Cauchy problem. Separately we will give some estimates 
of Z for large values of t, just to see the qualitative behavior of Z. 

Theorem 4.1. Suppose that (j, G C 2 [0, 1], fi(a) = a u fii(a), fii(ct) >p>0,0<a<l, u>0. 
Denote by e a small positive number. The function Z is infinitely differentiable for t ^ and 
x 7^ 0. The following estimates hold for < t < T. 
If n — 1, then 

\D™Z(t,x)\ ^Cr 1 ^, \x\ < 5,0 < m < 3. 

If n = 2, then 

\Z(t,x)\ < CrMoglarl -1 , \x\ < e, 
\D™Z(t,x)\ < Cr l \x\- m , \x\ < e, 1 < m < 3. 

Ifn>3, then 

\D™Z(t,x)\ < Cr 1 \x\- n+2 - m , \x\<e,0<m<3 

In all cases, 

\D™Z(t,x)\ < Ce- alxl (a > 0), \x\ > e _1 . 
The estimate ofD^Z, uniform in t, is as follows: 

| (D^Z) (f , x) | < C\x\- n - 2 e- a W (a > 0), |x| ^ 0. 

If \x\ < 5, £/ien 

|(D^Z) < Cr 2 |a:r ri+2 . (4.13') 



(4.8) 

(4.9) 
(4.10) 

(4.11) 
(4.12) 

(4.13) 



Proof. As before, using Jordan's lemma we write 
Z(t,x) = (27r)-f|x| 1 -f y e p */C(p)(p/C(p))5(i- 1 )K f _ 1 (|x|Vp^))^, x + 0. (4.14) 

The integral in (4.14) consists of the ones on T 7)W and . Let us begin with the first of them, 
denoted by Z°(t,x). Below we assume that 7 > e. 

If p G T 7jW , then p = •ye lip , \<p\ < un, \ < uj < 1. Under our assumptions, 

1 

*M = /^°"*(«)*.. 
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Let us consider the location of values of ptC(p), p G T 7jW . If \<p\ < ir/2, then Rep/C(p) > 0. 

1 

Suppose that ~ < tp < un. Then Rep/C(p) > _Rcos(u;7r), R = f 7 Q a !/ /i 1 (a) da, 

o 

1 </? vr/2 

Im p£(p) > p J a v sin(a^) da = ppr x ~ v J f3 v sin /3 dp > p(cutt)- 1 ^ J (3 V sin (3 d(3 > 0, 



so that < argp/C(p) < 7r, as p belongs to the part of T 7j£tJ lying in the upper half-plane. 
Similarly, — tt < argp/C(p) < for the part from the lower half-plane. Thus, | argp/C(p)| < 7T, 
and since T J>U is compact, we have | argp/C(p)| < <p < 7r, p G T, 



This means that 



Re \J p1C(p) > cos ^ ■ inf 



p a fj,(a) da 



dcf 



r > 



l l 
because Im J p a p(a) da = with p e T 7)(J only if p = 7, and there Re f 7 Q /i(a) <ia 7^ 0. 


Therefore, using the above-mentioned asymptotics of the McDonald function, we find that 



\Z°(t,x)\ < Ce-^l (a > 0), \x\ > e~\ 



As \x\ < e, we get 



C, 



if n = 1; 



|^°(M)| < S Clog|x|-\ if 71 = 2; 

(7|x|-" +2 , ifn>3. 



(4.15) 



(4.16) 



Let Z ± (t,x) be the parts of Z(t,x) corresponding to the integration over . If, for 



example, n > 3, then 

\Z ± (t,x)\ < C\x\ 



l-n 
1 



tr cos(oj7r) 



1 



log r \ log r 



n — 3 

/ \ 4 e -a[«[( 



_r_W2 

logrj a > 0. 



(4.17) 



1/2 



Let us make the change of variables z = z(r) = ( ) . In order to express (asymptot- 

\\ogrJ 

ically) r as a function of z, we denote s = logr, so that s~ 1 e s = z 2 where s — > 00 and z — > 00. 
Taking the logarithm of both parts of the last equality we get s — logs = 2 log z. It is known 
([32], page 50) that 



(4.18) 



s = 2 log z + O (log log z), z — > 00. 
Therefore r = r(^) satisfies the inequalities 

2 2 (logz)- fe < r(z) < (\ogz) b 

for some 6 > 0. 
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For |x| > the factor e * rcos ^) in (4.17) can be estimated by 1, and after the use of (4.18) 
the power terms, as well as the logarithmic ones, are dominated by the exponential factor (the 
integral in z is taken over (71, 00), > 0), so that 

\Z ± {t,x)\ < Ce~ a ' lxl , a'>0, 

and, together with (4.15), this implies (4.12) for n > 3, m = 0. 

For \x\ < e, the factor e a|x| uo g rJ j s estimated by 1, and an elementary estimate gives 
that 

\z ± (t,x)\ < cr 1 ^-"* 2 , 

which implies the required estimate of Z. 

The bounds for the derivatives, as well as the estimates (4.8)-(4.10) for n — 1 and n = 2 
are obtained in a similar way. Some of the estimates can in fact be slightly refined (using the 

t- 1 



asymptotic formulas for the Laplace integrals with logarithmic factors |32|), involving 

logt -1 

for small values of t, instead of t~ x . 

Let us prove (4.13). Let n > 3; the cases n = 2 and n = 1 are similar. The estimates of the 
McDonald function for small and large arguments can be combined as follows: 

|i^ f _x(^)| < C|^r^ +1 e-" 1 " 1 , z^O, (4.19) 

with a possibly different choice of the constant a > 0. 

Next, let us write an integral representation of O^Z. If u(t) = e pt , then 

t t 
(B M u) (t)=p J k(t - r)e pr dr = pe pt J k(s)e- ps ds. 


Using the expression (2.4) for k(s), the identity 

s~ a e- ps ds = p a - 1 -f(l-a,pt) 

([28]. 1.3.2.3), where 7 is an incomplete gamma function, we find that 

1 

(D^u) (t) = e pt J P a r ff a) 7(l - oc.pt) da. 


It is known that 

1 (l-a,z)~-^z 1 - a , z^O, 
1 — a 

7(1 — a, z) ~ T(l — a) — z~ a e~ z , z — » 00 
(see Chapter 9 in pQ). This implies the inequality 



7(1 ~ a,z) 
T(l-a) 



< C 
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valid, in particular, for all z = pt, p G S ltUI , t G (0, T]. Recalling also that 

p a fx(a) da = pfC{p) 



we see that the application of to the integral representing Z leads to the appearance of 
the factor \pJC(p)\ in the estimates of D^Z, compared to those of Z. Using also (4.19) and 
estimating by 1 the decaying exponential involving t (in the integrals over r~ w ) we come to the 
inequality (4.13). 

The proof of the inequality (4.13') is similar to those for estimates of the derivatives in 
spatial variables. ■ 

4.2. Subordination and positivity. Let us find a connection between Z and the fun- 
damental solution of the heat equation. Our approach follows [I] where the case n = 1 was 
considered (without a full rigor). 

Let us consider the function 

g(u,p) = JC(p)e- upK(p) , w>0,Rep>0. 
Let p = 7 + it, 7 > 0, r G 1R. As |r| — > oo, 

/C(p) ~ = , argp — > ±— . 

l°g V 7 2 + t2 + * argp 2 



It follows that 



Re(p£(p)) ~ C { ^== + 1- M 



log ^7 2 + r 2 2 (log a/t 2 + t 2 ) 



2 



whence 



e-^W I < C exp -an —J ■ - + — -J — , a > 0. (4.20) 



l og ( 7 2 + r 2) ( log ( 7 2 + r 2))2 

Writing log(7 2 + r 2 ) < C(7 2 + r 2 ) e , < s < 1/4, we find from (4.20) that 

oo oo 

J \g(u,j + iT)\dT<C J exp (-au [ j^^y + J^^ j) dr 





1 oo 



< C J e au W&? dr + C j e au rf+W° dr < Ce a " (7W + J e~ a ' 



dy 



<C + C I e~ a ' uzl 4e dz 
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(a' > 0), whence 

oo 

sup / \g{u, 7 + it) I dr < oo. (4.21) 

7>1 J 



-oo 



It follows from (4.21) (see [7]) that g(u,p) is the Laplace transform of some locally integrable 
function G(u,t): 

oo 

g(u,p)= [ e-*G{u,t)dt, (4.22) 



and the integral in (4.22) is absolutely convergent if Rep > 1. 

On the other hand, the function pJC(p) is positive and has a completely monotone derivative, 
so that e~ uplc ^ is completely monotone. Since fC(p) is completely monotone, we find that g 
is completely monotone in p (we have used Criteria 1 and 2 of the complete monotonicity; see 
[13]), so that G(u,t) > by Bernstein's theorem. 

Theorem 4.2. (i) The fundamental solution Z(t,x) satisfies the subordination identity 

oo 

Z(t, x) = J G(u, t)(47iu)- n/2 e- ]J £ du, x £ 0, t > 0, (4.23) 
o 

where G(u,t) > and 

G(u,t)du=l. (4.24) 



(ii) For all t > , x ^ ; Z(t, x) is non-negative, and 

J Z(t,x)dx = 1. (4.25) 



oo 



Proof. In order to prove (4.24), we integrate (4.22) in p using Fubini's theorem. We get 

oo 

dt / G(u, t)du= -, 
J P 



which implies (4.24). 

Let us prove (4.23). The convergence of the integral at infinity follows from (4.24), while 

near the origin the function u t— > (4iru)~ n > 2 e~~ decays exponentially. Let v(t,x) be the 
right-hand side of (4.23). Multiplying by e~ pt and integrating in t we find that 

oo 

pt v(t,x) dt = K{p) [ e- uplc{p \4iru)- n/2 e-^ du. 
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By the formula 2.3.16.1 from [28J, the right-hand side coincides with the one from (4.5), so 
that v(t,x) = Z{t,x). Now the non-negativity of Z is a consequence of (4.23), and the identity 
(4.25) follows from (4.23), (4.24) and Fubini's theorem. ■ 

4.3. Long time behavior. Let us give a rigorous proof of the asymptotics of the mean 
square displacement, basic for applications of the distributed order calculus. We also give some 
long time estimates of the fundamental solution Z. 



Theorem 4.3. (i) Let 

lffi(0) ^ 0, then 
If 

then 



m(t) — J \x\ 2 Z(t, x) dx. 

m{t) ~ Clogt, t — > oo. 
fi(a) ~ aa u , a — > 0, a, v > 0, 



m(t) ~ C(logt) 



t — >• OO. 



(4.26) 
(4.27) 
(4.28) 



(ii) Suppose that (4-27) holds with v > 1, if n = 1, and with an arbitrary v > 0, ifn>2. 
Then for \x\ < e, e > 0, and t > e^ 1 , 



C(logt)- — , i/n = l; 

Zit..r)< ^qbg|z||(log*)- , 'bg(logt), tfn = 2; 

> C|x|- n+2 (logt)- 1 '- 1 , ifn>3. 

Proof, (i) It follows from the Plancherel identity for the Fourier transform that 

m(f) = -(27r)" (A c Z(t,0) • 
Applying the Laplace transform in t we find that 



5=0 



(4.29) 



and after an easy calculation we get 



m{p) 



2n ■ (2tt) 
p 2 JC(p) 



whence 



m(t) = 2n • (27r)" y x(r) dr 
o 

where x was introduced in Sect. 3.1. Now the relations (4.26) and (4.28) are consequences of 
Karamata's Tauberian theorem IT3l. 
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(ii) As before, we proceed from the integral representation (4.14) where the contour 5 7)W 
consists of a finite part T 7iW and the rays r^ w . Let n = 1. Then (4.14) takes the form 

Z(t,x) = - [ e^-^Le-WV^Wrfp. (4.30) 

S7lj 

Asp — > 0, \/pfC(p) ~ CQogp" 1 ) - ", — ~ Cp _1 (logp~ 1 ) _ ^ ii , where ^-jp > 1. These 

vp£(p) 

asymptotic relations make it possible to pass to the limit in (4.30), as 7 — > 0, substantiating 
simultaneously the convergence to of the integral over T 7)£J and the existence of the integrals 
over the rays starting at the origin. 
Thus, 

00 

Z{t,x) < C J e rt««(^) r -i| logrl-^e-"^ 1110 ^ 1 " 1 ^ dr. (4.31) 
'0 

Let us decompose the integral in (4.31) into the sum of the integrals over (0, 1/2) and (1/2, 00). 
Estimating the latter we drop the factor containing \x\ and obtain easily the exponential decay, 
as t — > 00. The integral over (0, 1/2) is estimated via the function 



1/2 



1 



2 



M(t) = I e - art r- x ( log 3 J dr. 
Integrating by parts we see that 







/ 1/2 > 

M (t) < C I e"T + t J e - art ^log " rfr 



It is known (see (18.52) in [33J or (32.11) in [M]) that 

1/2 i-j, 

e- ari ^log^ 2 cir < Cr 1 (log t)^r 

for large values of t. This implies the first inequality of (4.29). 
Let n = 2. Then 

Z(p,x) = ^K(p)K (\x\y/jK&j, 
Ztt 

so that we have, for \x\ < e and small \p\, that 

Z(p, x) < C|p| _1 | log |p|r 1_I/ (log \x\~ l + log log IpT 1 ) . 

This estimate is sufficient (for v > 0) to substantiate passing to the limit, as 7 — * 0. The 
above argument gives, as the main part of the upper estimate of Z(t,x) for a large t, the 
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expression 

1/2 



C^oglx]- 1 ) J e- art r- 1 \\ogr\- l - v (loglogr-- 1 ) dr 



o 

1/2 



d (log^r 1 ) J e~ art (loglogr- 1 ) — (logr" 1 ) " dr. 



Integrating by parts we reduce the investigation of the above integral in r to that of two 
integrals, 

1/2 



J e~ art r~ l (logr- 1 ) 1 " dr 



o 



(it has been estimated above), and 

1/2 



J e~ art (logr" 1 ) "loglogr -1 ^ ~ Cr 1 (log t)-"log(log*), t -> oo 
o 

3^] . 32.11). This results in the second estimate from (4.29). The third one is derived similarly. 



The relations (4.26) and (4.28) for the case where n — 1 and fj,(a) = const or /x(a) = const-a^ 
were proved in [4]. 

5 The Cauchy Problem 

5.1. The homogeneous equation. Let us consider the equation (4.1) with the initial 
condition 

u(0,x) = <p(x), xeW 1 , (5.1) 
where (p is a locally Holder continuous function of the sub-exponential growth: for any b > 0, 

\<p(x)\ <<V N - (5.2) 

We will assume that the weight function /x defining the distributed order derivative D^' satisfies 
the conditions of Theorem 4.1. 

Theorem 5.1. (i) The function 

u(t,x) = J Z(t,x-0<p(0d£ (5.3) 
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is a classical solution of the Cauchy problem (4-l)-(5.1), that is the function (5.3) is twice 
continuously differentiable in x for each t > 0, for each x E M. n it is continuous in t > 0, the 
function 



1 1— > J k(t — t)u(t,x) dr, t > 0, 



is continuously differentiable, the equation (4-1) is satisfied, and 

u(t,x) — > (p(x), as t — > 0, (5.4) 

for allx G W 1 . 

(ii) On each finite time interval (0,T] 7 the solution u(t,x) satisfies the inequality 

\u(t,x)\ < Ce dlxl , i6l", (5.5) 
with some constants C, d > 0. If tp is bounded, then 

\u(t,x)\<C, xeR n , 0<t<T. (5.6) 

(iii) For each x EMJ 1 , there exists such an e > that 

\B M u{t, x) | < C x t~ 1+£ , < t < T. (5.7) 

Proof. Using (4.25) we can write 

u(t,x) = J Z(t,x-S)[<p(Z)-<p(x)]dS + <p(x). (5.8) 

Let us fix x and prove (5.4), that is prove that the integral in (5.8) (denoted by uo(t,x)) tends 
to 0. 

Let n = 1. Then 

7+ioo 

1 f „vt K(P) 



7— 100 

where 7 > 0, 



Mt - x)= ^J. ert vm H(p ' x)dp (5 - 9) 



(the change of the order of integration leading to (5.9) will be justified when we prove the decay 
of H (p, x), as p — > 7 ± ioo; see below). 
By our assumption, 

\<p(x) - <p(£)\ <C x \x-S\\ A>0, \x-£\<l. 
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Let p = \fpiC{v)i p = 7 + it. As \t\ — > oo, p ~ C I - — -— I . We have 

V.log|r|y 

\H(p,x)\<C J e- plx -t\\x-Z\ x dZ + C J e - plx -^ +b ^ d£ + \(p(x)\ J 

\x-£\<l \x-£\>l \x-£\>l 

oo 



< Cp- X - X + Ce blxl J e (b - p)lzl dz + 2\<p{x)\ J e~ pz dz < Cp- 1 -\ 



\z\>l 1 

if b is taken such that p > b. Therefore the absolute value of the integrand in (5.9) does not 
exceed 

Ce^lrr'-^loglrl)^ 2 , 
so that the integral in (5.9) exists and possesses a limit, as t — > 0, equal to 

7+ioo 



7— lOO 



The integrand in (5.10) is analytic in p on the half-plane Rep > 7. Let us consider (within 
that half-plane) a contour consisting of an interval {p : Rep = 7, \p\ < R} and the arc 
{p : Rep > 7, \p\ = R}, R > 7. The absolute value of the integral over the arc (with the same 
integrand as in (5.10)) does not exceed CR~ x / 2 ( log R) x ^ 2 — >■ 0, as R — > 00. This means that 
m (0, x) = 0, and we have proved (5.4) for n = 1. The scheme of proof is completely similar for 
n > 1 too; one has only to use the asymptotics of the McDonald function. 

If we perform the above estimates, not ignoring the dependence on x but, on the contrary, 
taking it into account, then we obtain the estimates (5.5) and (5.6). 

Due to the estimates of Z given in Theorem 4.1, we may differentiate once in the spatial 
variables in (5.3) under the sign of integral. Using also the identity (4.25) we get, for a fixed 
x°, the formula 

Let us decompose the domain of integration in (5.11) into the union of 

^ = {(£8": \x° -C| > 1} 

and fi 2 = IR n \ Correspondingly, the integral becomes a sum of two functions, Wi(t,x) + 
W2(t,x). If x is in a small neighbourhood of x°, while £ G Q\, then \x — £| is separated from 
zero. Therefore 

d Wl (t,x ) f d 2 Z(t,x°-Q f o 

— ^— = y — ^ — b(o - )] de. (5.12) 
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Let d be a small positive number, d — (0, . . . , d, . . .), with d being at the k-th place. Then 
I - „ , „,- fd 2 Z(t : x° O 

n 2 



Ox 2 



-M0-<p(x°)]dt; 



w 2 (t,x° + d) -w 2 (t,x°) - J 

I 



\x°-£\<2d 



\x°-£\<2d 



+ 



/ w 



2cK|:r -£|<l 



\x°-(]<2d 

8Z(t,x° + az(t,x°-o 



dx k 



dxi 



d 2 Z(t,x°-0 



dx 2 



[tp(0 - v(x )] dC 
(5.13) 



The integrals converge due to the local Holder continuity of (p. 
We have (if n > 2) 



\ 1 



dZ{t,x° + d-£) 
dx k 



wo - v^ )] <*e < ct-'d- 1 J \x° + d - er n+i ie - *°i A « 



|z°-£|<2<i 



|z°-£|<2d 



= Cf 



"M -1 J \r] + d\- n+1 \r]\ x dr] < Cr 1 d x -> 0, d^O 



r|<2d 



(the change of variables r] = d( was made in the last integral). In a similar way we obtain 
estimates of other integrals over the set {\x° — £| < 2d}. 

In the integral over its complement, we use the Taylor formula: 



1 

d 



dZ(t,x° + d-g) dZ{t,x°-£) 



dx k 

If \x° - f | > 2d, then 



dxk 



d 2 Z{t,x°-£) dd 3 Z{t,x° + 0d-f) 



, < < 1. 



1 



Using the inequality for the third derivative of Z from Theorem 4.1 we find that the last integral 
in (5.13) does not exceed 



Cdt 



- 1 J \i - x Q \- n ' l+x di < Ct~ 1 d x -> 0, 



2d<|x°-£|<l 



as d — > 0. 

It follows from (5.12), (5.13) and the above estimates that 



d 2 u(t,x°) fd 2 Z(t,x°-0 



dx\ 



I 



dx k 



[<p(Z) - V (x )} dc 



(5.14) 
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If n = 1, then the formula (5.14) is obtained by a straightforward differentiation under the sign 
of integral. 

Let us consider the distributed order derivative O^w. First of all we check the identity 

B M Z(t, x) = AZ(t, x), t > 0, x ^ 0. (5.15) 

A direct calculation based on identities for the derivatives of the McDonald function [1] shows 
that AZ(p,x) = pJC{p)Z{p,x). On the other hand, if x ^ 0, then Z(t,x) — > 0, as t — > 0. 
This fact follows from the integral representation of Z in a manner similar to the above proof 
of (5.4). Therefore the Laplace transform of ~B^Z(t, x), i / 0, equals p)C(p)Z(p,x), which 
implies (5.15). 

Now, having the estimates of the derivatives of Z in spatial variables given in Theorem 4.1, 
from (5.15) we get estimates for H^Z sufficient to justify the distributed differentiation in 
(5.8). Thus we come to the formula 

(B^u) (t, x°) = [ (DWZ) (t, x° - OHO ~ <p(x°)} dt (5.16) 



Together with (5.14) and (5.15), this proves that u(t,x) is a solution of the equation (4.1). 

In order to prove (5.7), we use the inequalities (4.13), (4.13'), and the assumption (5.2) with 
b < a. Substituting into (5.16) we get, for a fixed x°, that 



(DWu) (t, x°) | < Cr 2 J \x° - £\- n+2+x d£ + C J \x°- £|-«-2+A e -a|*°-£l ^ 

[xO-^|<tV2 tV2<[a;0_ ? [<i 

+ C J \x°- £p- 2 e -«l*°-*l (e fo l«l + e^°l) d£ 



\x°-£\>l 

jl/2 



< Cr 2 I r 1+x dr + C I r - 3+x e- ar dr + C < Cr 1+ ^ 



t l/3 

for small values of t, as desired. ■ 

5.2. The inhomogeneous equation. Let us consider the Cauchy problem 

(d^w) (t,x) - Au{t,x) = f{t,x); xeR n ,t>0, (5.17) 

u(0,x)=0. (5.18) 

We assume that the function / is continuous in t, bounded and locally Holder continuous in x, 
uniformly with respect to t. Our task in this section is to obtain a solution of (5.17)-(5.18) in 
the form of a "heat potential" 

t 

u(t, x)= I dr I E{t-r,x- y)f{r, y) dy. (5.19) 
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In contrast to the classical theory of parabolic equations |14j . the kernel E in (5.19) does not 
coincide with the fundamental solution Z - just as this happens for fractional diffusion equations 
pTOj, [9] . However the behavior of the function E is very similar to that of Z. Applying formally 
the Laplace transform in t and the Fourier transform in x we find that 

E(p,0 



pK{p) + K| 2 

whence 

E(p,x) = (27r)-^|x| 1 -t( J9 /C(p))5(f- 1 )K f _ 1 (|s| v / p/C(p)), (5.20) 

which differs from (4.5) only by the absense of the factor JC(p) with a logarithmic behavior at 
infinity. Therefore the function E(t,x) obtained from (5.20) via contour integration, satisfies 
the same estimates (see Theorem 4.1) as the function Z, except the estimates for large values 
oft. 

The function E(t,x) is non-negative. Indeed, the function p i— »■ p v l 2 K v {a^fp\ a > 0, is the 
Laplace transform of the function 

t I — > P 4t 

(2ty+ 1 

(see [7j). This means that the above function in p is completely monotone. Since the function 
p i — ► pJC{p) is positive and has a completely monotone derivative, we find that E(p,x) is 
completely monotone in p, so that E(t,x) > 0, i ^ 0. 

The counterparts of the estimates (4.29) (proved just as in Theorem 4.3) are as follows. If 
(4.27) holds with v > 0, then for \x\ < e, e > 0, and t > e~ l 

( Cr 1 (log t) 1 ^, if n = l; 

E(t,x) < I Cr 1 log log Hog ifn = 2; (5.21) 

(Cr l \x\- n+2 , ifn>3. 

The function E has (in x) an exponential decay at infinity. 

In fact, for the analysis of the potential (5.19) we need estimates of E and its derivatives, 
uniform in t e (0, T]. 

Proposition 5.2. Let // satisfy the conditions of Theorem 4-1- Then, uniformly in t G (0, T], 
\DiE(t,x)\ < C\x\- j ~ n \l + |log|x||| / V a|x| , x^0,j>0, (5.22) 

B^E{t,x) <C\x\- n ~ 2 \l + \log\x\\fe- alx \ x + 0, (5.23) 
whre C, a, (3 are positive constants. 

Proof. Let, for example, n > 3 (other cases are considered in a similar way). As usual, 
we write the Laplace inversion formula and deform the contour of integration to S' 7iU . The 
integral over T 7AJ gives an exponentially decaying contribution without local singularities. In 
the integrals over the rays we use the upper bound 

\Kn !(*)! <C\z\S +1 e- a W, z^0, 
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(a > 0) obtained from the asymptotics of the McDonald function near the origin and infinity. 



As in the proof of Theorem 4.1, we perform the change of variable z = and use the 

\logr / 

inequality (4.18) for the dependence of r on z. 

As a result, for the integrals over we obtain the upper bound 

oo 

C W -/,(lo g ^e-.|.l^,<C|,|-»(UogWI + lA-W 



71 

with some positive constants, and we come to the estimate (5.22), j = 0. The estimates of the 
derivatives in spatial variables are proved similarly. 

The proof of (5.23) is completely analogous to that of the inequality (4.13) for 3^Z. ■ 

As we have noticed, 

E(p,x) = ^-yZ(p,x), 

and since 

1 

Z(v, x) dx = -, 
V 



I 



we have 

so that we come to an interesting identity 

J E(t,x)dx = x(t). (5.24) 

The existence of the integral in (5.24) follows from the above estimates or from the fact that 
x G L l °° (see (3.2)) and Fubini's theorem. 

Theorem 5.3. Under the above assumptions regarding f , and the assumptions of Theorem 4-1 
regarding fi, the function (5.19) is a classical solution of the problem (5.17)-(5.18), bounded 
near the origin in t for each x Gl". 

Proof. The initial condition (5.18) is evidently satisfied. Just as for the kernel Z above, we 
prove that — AE = for x ^ 0. Next, we may differentiate once in (5.19) under the sign 
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if integral. Indeed (here and below we make estimates for n > 3; other cases are similar), 

t t 
dE(t — t,x — y) 



hi 



dxi 



iy <cf d r j W -( 1 + |.„ gW |)V»N d , 



\x\>^F 
t 

+ C f r^dr 

\x\<V* 



x\- n+1 dx 



t I 

<cJr- l l 2 dT J \y\- n - l (l + \og\y\ + ]^\ogT- l )dy + C J r- l l 2 dr J \y\~ n+1 dy < oo. 



Ivl>i 



■Ki 



In order to calculate the second order derivatives, note that the function 



t-h 



u h (t, x) = J dr J E(t - t,x - y)f(r, y)dy, t> h> 0, 



R" 

may be differentiated twice, and that 

d 2 



I 



dx\ 



E(t - r, x - y) dy = 0, 



whence 



d 2 u h (t,x) 
dx 2 



t-h 



J dr J JL E (t-T,x-y)[f(T,y)-f(T,x)} 



dy. 



(5.25) 



Using the local Holder continuity and boundedness of /, we perform estimates as above and 
prove the possibility to pass to the limit in (5.25), as h — > 0, so that 

t 

Au(t,x) = j dr J AE(t - r,x - y)[f(r,y) - f(r,x)]dy. (5.26) 

I" 

To calculate u, we use (5.24) and write 
u(t,x) = j 'It J E(t-T,x-y)lf{T,y)-f{T,x)]d v + j x(t-T)f(T,x)d T = Ui(t,x)+u 2 (t,x). 

OR™ 

Recall that u 2 {t,x) = (*, 

x), so that D^ti2 = / (see Sect. 3). 
dE 

Let us consider u\. First we estimate — — . As before, we use the contour integral represen- 
ts 

tation of E and note that the differentiation in t leads to an additional factor p in the integrals. 
This results in the estimates 



dE(t,x) 



< Cr 2 \x\- n+2 , \x\ < e; 



(5.27) 
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dE(t,x) 



dt 



< C , |x|^ 2 (|log|x|| + l) /3 e- a|:r| , \x\^0. (5.28) 



As the first step of computing D^'ui, we compute — — . Note that 

at 



7+ioo 



where 



J E (t -r,x- y)[f(r, y) - f(r, x)] dy = ^y /2+l ■ J e^\pK{p))^-^ L n (p,x,r) dp 

7-ioo 

(5.29) 

L n (p, x,t)= J \x- y] 1 '^ [f(r, y) - f(r, x)\K^{\x - y\\/pK,{p)) dy. 



The role of the function L n is quite similar to that of the function H introduced in the proof 
of Theorem 5.1 (where the case n = 1 was considered in detail). Using, as it was done there, 
the local Holder continuity and boundedness of / we find that 

\L n {p,x,r)\<C\^plC{p)\^- x -\ 

As in the proof of Theorem 5.1, we deform the contour of integration to the right of the line in 
(5.29) and show that 

Jim J E(t - T ,x - y)[f(r,y) - f(r,x)]dy = 0. (5.30) 

R" 

On the other hand, using (5.27) and (5.28) we get 

/ 9E{t ~£ X " V) |/(r, y) - /(r, x)\ dy < C J \y\~ n - 2+A (| log |y|| + if e-M dy 



\y\>Vt^ 

+x dr 



+ C(t-r)- 2 J \y\- n+2+x dy = 2C J r- 3+A (| logr| + lfe~ ar dr + 2C(t - r)~ 2 J r 1 

< C(t - r)- 1+A / 2 (| \og(t -t)\ + If. 

Together with (5.30), this implies the equality 

t 

dui f f dE(t — r,x — y) 



dt 



JdrJ dE ^^- y \ f { r,y)-f(r,x)]dy. 



o 



Now we compute H^ui using the formula (2.3), the fact that k G L l ° c (following from (2.4)) 
and Fubini's theorem: 

t 

(D M «i) (t,x) = J dr J (BM E ) (t-T,x-y)[f(T,y)-f(r,x)}dy. 

R n 
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Together with (5.26), this means that Am = D^iii = D^w — /, as desired. ■ 

In Theorem 5.3 we constructed a solution u of the problem (5.17)-(5.18), such that u = 
ui + ii2, ui(0,x) = 112(0, x) = 0, Mi is absolutely continuous in t, and U2 = I^'f. On this 
solution u, 

I M D W « = I&>(k * u[) + I (At) / = x * k * u[ + u 2 = ux + u 2 = u 
(u' means the derivative in t). Applying iw to both sides of the equation (5.17) we find that 

u(t,x) — J x(t — s)Au(s, x) ds = (x * f)(t, x). (5.31) 

Q 

The equation (5.31) can be interpreted as an abstract Volterra equation 

u + x * (Au) = ip, (5.32) 

if we assume that u belongs to some Banach space X (in the variable x), and A is the operator 
—A on X. The operator —A generates a contraction semigroup if, for example, X = L 2 (M. n ) 
or X = Coo(M n ) (the space of continuous functions decaying at infinity; see Sect. X.8 in 
[3 lj ) . Now the existence of a solution in Li(§,T\X) can be obtained from a general theory 
of equations (5.32) developed in [6]; it is essential that x is completely monotone (conditions 
of some other papers devoted to equations (5.32) do not cover our situation). Of course, our 
"classical" approach gives a much more detailed information about solutions, while the abstract 
method is applicable to more general equations. 



6 Uniqueness Theorems 

6.1. Bounded solutions. In this section we consider a more general equation 

(0 ( ^w) (t, x) = Lu{t, x), x e W 1 , < t < T, (6.1) 
with the zero initial condition 

u(0,x)=0. (6.2) 

Here L is an elliptic second order differential operator with bounded continuous real-valued 
coefficients: 

Lu=J2 Oij (t, X )-Q^r. + Yl h i (*> x ) -fa. + c (*> 

n 

J2 <kj{t, x)tej > 0> + i = ■ ■ ■ >U e M n . 

We assume that « G C 3 [0, 1], u(l) 7^ 0. 

We will consider classical solutions u(t, x), such that (B^'u) (t, x) belongs, for each fixed x, 
to L p (0, T) with some p > 1. As we saw in Theorem 5.1 and Theorem 5.3, the solutions for the 
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case L = A obtained via the fundamental solution and heat potential possess the last property 
making it possible to represent D^'u in the Marchaud form (3.9). 
It is often convenient to transform the equation (6.1) setting 

u(t,x) = u\(t)w(t,x) 

where A > 0, and u x is the solution of the equation ~B^u x = Xu x constructed in Sect. 2.3. It 
is easy to check that the function w satisfies the equation 

(A x w) (t,x) = (L - X)w(t,x) 

where 

(A x w) (t,x) = — - I k(t)w(t,x) +lim / ux{r)k\t-T)[w(T,x) - w(t,x)]dr > . (6.3) 
u\{t) J ' | 

The operator (6.3) is very similar in its properties to the distributed order derivative D^. 

Theorem 6.1. Ifu(t,x) is a bounded classical solution of the problem (6.1)-(6.2), such that 
for each x E R n , B^u E L p (0, T) for some p > 1, then u(t, x) = 0. 

Proof. Let M = sup \u(t,x)\. Consider the function 

F R (t, x) = ^\x\ 2 + a J x(s) ds + 1 j , 
t 

with R, a > 0. It follows from (3.2) that f x(s) <is — > 0, as £ — > 0. As we have seen (Sect. 3.2), 

o 

Q x(s)ds^ = 1, so that 

(DMfi«)«,») = ^. 

Let c = sup |c(£,x)|, d > 0, X = c + d. Since «a is non-decreasing (Theorem 2.3), and 
k'(s) < 0, we have 



u x (T) R 2 u x (TY 
On the other hand, 

2M ( n n \ 2M 

(LF R ) (t,x) < -jp ^a^(x) < -jp (Ci + C^l) , Ci,C 2 > 0, 

\i=l j=l / 

so that taking d>0we get 

((A x — (L — A)) F fl ) (£, x) > ^ - 2d - 2C 2 \x\ + rf|a;| 2 + d) > 

i? 2 W(T) y 
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for all x G IR n , t G (0, T), if a is taken sufficiently big. 

Denote t>(t,x) = u(t,x) — Fn(t,x). By the above inequalities, 

(A A f ) (t, x) — (L — A)f (t, x) < 0. (6.4) 

If |x| = R, then v(t, x) = u(t, x)-M-R~ 2 (^a f x(s) ds + l^j < 0. Next, v(0, x) = -F R (0, x) < 

for all x. This means that v(t,x) < for |x| < R, t G [0, T]. Indeed, otherwise the function 
v would possess a point of the global maximum (t°,x°) on the set {(t,x)\ < t < T, \x\ < R}, 
such that v(t°,x°) > 0. Then 

(L - X)v(t°,x°) < 



(see the proof of the maximum principle for a second order parabolic differential equation [111 
ED]), so that ( A xv) (f,x°) < 0, due to (6.4). However it follows from (6.3) that (A\v) (t°,x°) > 
0, and we have come to a contradiction. 
Thus, we have proved that 

u(t, x) < — ( \x\ 2 + a I x(s) ds + 1 | , \x\ < T. 

Since R is arbitrary, we find that u(t,x) < for all t G [0, T], x G M n . Considering —u(t,x) 
instead of w(t, x) we prove that u(t, x) = 0. ■ 

The above proof was based on standard "maximum principle" arguments. In fact, it is 
easy to prove, for the equation (6.1), an analog of the maximum principle itself. Namely, let 
c(t, x) — A < for t G [0, T], x G G, where G C M n is a bounded domain. Suppose that 

(L - X)u(t, x) - {A x u) (t, x) > 0, (t, x) G [0, T] x G. 

Then, if sup u > 0, then 

[0,T]xG 

sup u = sup u. 

[0,T]xG [0,r]x9G 

The proof is similar to the classical one [20] . 

6.2. Solutions of subexponential growth. In this section we will prove a more exact 

d 2 

uniqueness theorem for the case where n = 1, L = — — 

ox 1 

Theorem 6.2. Suppose thatu(t,x) is a classical solution of the problem (6.1)- (6. 2) with n = 1, 
d 2 

L = — — such that for any a > 0, 
ox 2 

\u(t,x)\ < C a e aN , 0<t<T, xGM 1 , 
and D, u G L p (0, T), p > 1, in t, for any fixed x. Then u(t, x) = 0. 
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Proof. This time we choose the comparison function as 

F£\t,x) = Me aR [Z(t,x - R) + Z{t,x + R)}, \x\ < R, 



(6.5) 



where Z is the above fundamental solution of the Cauchy problem (Sect. 4), M and a are 
positive constants to be specified later. We will need the following auxiliary result. 

Lemma 6.3. For any T > 0, there exists such a constant p > that 

Z(t,0)>p o , 0<t<T. 



Proof. By (4.7), 



z{ P ^) = l -^ P ~ z mi) 1 



so that 



Z(p,0)^^M{plogp)-^, 



p — > oo. 



(6.6) 



Note that in Sect. 4 we used the Laplace inversion formula for Z(t, x) only for i / 0. Here the 
task is just the opposite, and we use the inversion formula from [11] involving the derivative of 
the Laplace image. In our case 



dp 



so that 



fc'(p) = J ( a ~ l)P a 2 A i ( a ) da = o{p 1 ), p — » oo, 
o 

dZ(p,0) 



dp 



< c e \p[ 



Rep > 1, 



for any e > 0. This is sufficient for the inversion formula 

7+ioo 

Z(t,0) = ^- J Z(p,0)eP t dp, t^O, 

7— ioo 



(6.7) 



where 7 > 1. Using (6.6), (6.7) and an asymptotic theorem for the Laplace inversion (see 
(22.115) and (22.114) in [33]) we find the asymptotics 



Z(t,0) = Cr 1 / 2 (logi) 



-3/2 



t -> +0. 



(6.8) 



For our purpose, it is sufficient to derive from (6.8) that Z(t, 0) — > +00, as t — > +0. 

On the other hand, it follows from the subordination identity (4.23) and the fact that 
G(u,t) > for each t on the set of a positive measure in u (see (4.24)), that Z(t, 0) > for 
each t > 0. Together with (6.8), this implies the required inequality. ■ 
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Proof of Theorem 6.2 (continued). If \x\ = R, that is x = ±R, then by (6.5) and Lemma 

6.3, 

F^(t,x) = Me aR [Z(t,0) + Z(t,±2R)] > Mp e aR . 

We have u(t, x) < F^\t, x), if C a < Mp (a has not yet been chosen), that is if M is chosen in 
such a way that M > C a p^ 1 . 

The function w(t,x) = F^\t,x) — u(t,x) satisfies, if \x\ < R, t G (0, T), the equation 

B^w = If \x\ = R, then w(t,x) > 0. In addition, 



dx 2 



liminf w(t, x) > 0, 



if \x\ < R. It follows that w(t,x) > for t G (0,T], |x| < i?. 

Indeed, if w(t,x) < for some t and x, then there exist t° G (0,T], a; G M 1 , < -R, such 
that 

w(t°,x°) = inf w(t,x) < 0. 

|x|<R 

te(o,T] 

If \x\ < R, then the function is infinitely different iable in t, with the derivative being 
continuous on [0, T\. Therefore we may write W^w in the Marchaud form, so that 

t 

k(t)w(t, x) + lim [ fc'(r) \w(t -t,x)- w(t, x)} dr = - ^f^ - (6-9) 

e^0 J dx 2 

e 

For (t,x) = (t°,x°), we see that the left-hand side of (6.9) is negative, while the right-hand 
side is non-negative, and we get a contradiction. Thus, we have proved that 

u(t,x) <Me aR [Z(t,x-R) + Z(t,x + R)}, < t < T, \x\ < R. (6.10) 

Let us fix x and consider the limit R — > oo. For large values of R we have 

Z(t,x±R) < Be- bR , tG(0,T], 

where b > depends only on T, B > depends on T and x, and not on i?. By (6.10), 

u(t,x) < BMe ia - b)R . (6.11) 

Now we choose a in such a way that a < b, and (see above) fix M > C a p 1 . Obviously, M 
does not depend on i?. Passing to the limit in (6.11), as i? — > oo, we see that u(t,x) < for 
arbitrary t and x. Similarly, taking —u(t,x) instead of u(t,x), we find that u(t,x) > 0. ■ 
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